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ABSTRACT 

It is often suggested that the distant galaxies recently identified in 850-^tm surveys 
with the SCUBA bolometer array on the JCMT Telescope are high-redshift analogues 
to local ultraluminous infrared galaxies, based on their similar spectral energy distri- 
butions and luminosities. We show that these two populations of objects must differ 
in at least one fundamental way from each other. This assertion is based on a con- 
sideration of the possible fates of gas in the high redshift SCUBA galaxies, given the 
requirement that they most evolve into some subset of the low-redshift galaxy pop- 
ulation with a comoving density of about 1CP 4 Mpc~ 3 . One possibility is that the 
SCUBA galaxies have similar gas density profiles to local ultraluminous galaxies. If 
this is the case, then they must derive almost all their power from AGNs, which ap- 
pears not to be the case for local ultraluminous galaxies, which are predominantly 
star-formation powered. Another possibility is that the SCUBA galaxies have more 
extended gas density profiles than local ultraluminous galaxies. In this case they must 
be almost all star-formation powered, and much of the star formation in the Universe 
can happen in these objects. Either way there is a significant difference between the 
low- and high- redshift populations. 

Key words: galaxies: formation - infrared: galaxies - quasars: general - cosmology: 
observations 



1 ULTRALUMINOUS GALAXIES AT LOW 
REDSHIFT 

Ultraluminous infrared galaxies (ULIGs) are peculiar galax- 
ies with infrared luminosities in excess of 10 12 Lq (Sanders 
& Mirabel 1996). These galaxies have very dense gas cores, 
as indicated by molecular line measurements, both of CO 
(Sanders, Scoville & Soifer 1991), and also of high-density 
traces like HCN (Solomon, Downes & Radford 1992) and CS 
(Solomon, Radford & Downes 1990). A number of measure- 
ments at radio, infrared, and optical wavelengths (Condon 
et al. 1991, Surace et al. 1998, Surace & Sanders 1999, Soifer 
et al. 1999) suggest that most of the very substantial bolo- 
metric luminosities of these galaxies comes from these dense 
molecular cores. 

The existence of these extreme objects leads to the ques- 
tion: what do they evolve into? One approach to answering 
this question is to look for stellar systems at redshift z = 
with central densities (see Fig. 1 of Binggeli 1994 for measre- 
ments of stellar densities in local galaxies) of about 100 Mq 
pc -3 or greater, the gas densities of the ULIG cores. Kor- 
mendy & Sanders (1992) made this comparison and showed 
that the only realistic candidates on these grounds at z = 



are giant elliptical (gE) cores. Galaxy disks (which have 
lower surface-brightnesses than gE cores; Freeman 1970) do 
not have high enough stellar densities. These gE cores are 
dense but have profiles that are flatter than an extrapolation 
of a de-Vaucouleurs r* law to small radii; they are rarely 
isothermal (Lauer 1985). Recent Hubble Space Telescope 
imaging (Faber et al. 1997) has confirmed the existence of 
such cores in the most luminous ellipticals (they are, how- 
ever, less common in lower-luminosity ellipticals). Often the 
cores are kinematically decoupled from the rest of the galaxy 
(e.g. Forbes et al. 1996) , suggesting that their formation hap- 
pens in a different way to that of the rest of the galaxy. These 
cores also tend to have supermassive black holes (BHs) in 
their centers, as inferred from stellar-kinematical measure- 
ments. The bigger the galaxy, the bigger this BH (Magorrian 
et al. 1998, van der Marel 1999, Ferrarese & Merritt 2000, 
Gerbhardt et al. 2000). Additionally, the bigger the galaxy, 
the bigger the core (Lauer 1985). 

An additional connection between the dense molecular 
cores of the ULIGs and the gE stellar cores + BHs is that 
their masses and sizes are similar, in addition to their densi- 
ties. Typical masses for both are several times 10 9 Mq and 
typical radii are several hundred parsecs (see Sakomoto et 
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al. 1999 for details of Arp 220, the nearest ULIG, and Lauer 
et al. 1985 for the core properties of gE cores). The BH con- 
tribution to the gE core masses can be added in from the 
correlations of Magorrian et al. 1998 or van der Marel 1999 
which relate the BH mass to the galaxy mass, in combina- 
tion with Tables 2 and 3 of Lauer 1985, which relate the 
core size to the size of the whole galaxy. Masses of the stel- 
lar cores can be computed directly from the core sizes, along 
with the mass-to-light ratio of Fukugita, Hogan & Peebles 
(1998). 

But locally, gE galaxies are much more common than 
ULIGs. The number density of ULIGs with far-infrared lu- 
minosities in excess of 10 12 L© is about 1.1 x 10 -7 Mpc -3 
(Saunders et al. 1990). Typical ULIGs have total Tf-band 
(2.2 /im) absolute magnitudes of Mk = —25.4 (Carico et 
al. 1988, 1990; here, as throughout this paper we assume an 
Einstein-de Sitter cosmology with Ho = 50 km s _1 Mpc -1 ), 
most of which comes from pre-existing stars in the pro- 
genitor galaxies (the dense cores are optically thick at K- 
band; Goldader et al. 1995). The number density of normal 
galaxies with M K < -25.4 is 1.6 x 10~ 4 Mpc~ 3 (Szokoly et 
al. 1998), a factor of about 1500 larger than the density of 
ULIGs. Most of these A~-band-luminous normal galaxies are 
early-type galaxies, which have the cores + BHs discussed 
in the previous two paragraphs; the contribution of lumi- 
nous disks, which tend to be blue, to the if-band luminosity 
function at the very bright end is small (see e.g. Binggeli, 
Sandage & Tammann 1988 for the optical luminosity func- 
tion of the different Hubble types, and Huang et al. 1997 for 
the optical-near- infrared colors) . 

Therefore if gE cores were made by ULIGs, the comov- 
ing density of ULIGs must have been a factor 1500 higher 
in the past than it is today, approximately 10~ 4 Mpc~ 3 . 
But this is very similar to the comoving density (e.g. Tren- 
tham, Blain & Goldader 1999) of the SCUBA (> 1 mjy 
at 850 /im in the observer frame) sources (Smail, Ivison & 
Blain 1997, Barger et al. 1998, Hughes et al. 1998, Eales et 
al. 1999), which are presumably high-redshift objects, and 
which Barger et al. (1998) argue are ULIGs based on their 
similar bolometric luminosities and spectral energy distribu- 
tion (SEDs). Therefore a picture in which SCUBA sources 
make gE cores and their associated black holes seems very 
likely not just based on similarities of their SEDs and lumi- 
nosities to those of local ULIGs but also due to the comoving 
number density of local gE cores + BHs being similar to that 
of the high-redshift SCUBA galaxies. It is this possibility 
that we investigate in the rest of this paper. 



2 ULTRALUMINOUS GALAXIES AT HIGH 
REDSHIFT 

In the previous section we argued that the SCUBA sources 
are probably the high-redshift analogues of low-redshift 
ULIGs, both in terms of their SED and luminosity concor- 
dance, and also in terms of the comoving number density 
concordance between the SCUBA galaxies and the likeliest 
remnants of the local ULIGs (the gE cores + BHs). 

This SED (Ivison et al. 1998) and luminosity (Barger 
et al. 1999) concordance is observed for individual sources, 
and is also implied on statistical grounds. By this it is meant 
that: if the SCUBA sources were significantly hotter than 



local ULIGs, then they would overproduce the far-infrared 
background shortward of about 500 /im measured by Fixsen 
et al. (1998) using COBE. If they were signifiantly colder 
than local ULIGs, they would severely underpredict the 
background at 450 /im: the SCUBA sources are observed to 
generate at least one-third of the background measured by 
Fixsen et al. at 450 /im (Blain et al. 1999a). These two state- 
ments assume a median source redshift of approximately 
3 (see the radio measurements of Smail et al. 2000) and 
then follow directly from the result that almost the entire 
background at 850 /im is generated by the SCUBA sources 
(Blain et al. 1999b, Barger, Cowie & Sanders 1999). Both 
statements are independent of any specific model of SCUBA 
source luminosity and density evolution. 

Now let us make the following additional assumption, 
which we will call Al: 

The density structure of the gas in the SCUBA 
sources is the same as in local ULIGs. 

This is an assumption, not an observation, since we 
have no direct probes of the gas densities in the SCUBA 
sources. We only know about CO in two of the SCUBA sor- 
ces (Frayer et al. 1998, 1999), and have no measurements 
at all of high-density tracers like HCN or CS in these ob- 
jects. There are two indirect pieces of evidence support- 
ing Al. Firstly, HCN has been detected in the Clover leaf 
quasar (Barvainis et al. 1997), a strongly-lensed infrared 
quasar which has an SED at far-infrared wavelengths sim- 
ilar to those of the SCUBA sources. Secondly, if the burst 
producing the high far-infrared luminosity is shortlived (see 
e.g. Blain et al. 1999a), then it is unlikely that it can be 
maintained in a coherent way over a galactic scale (in this 
case, the large scale over which the luminosity emerges from 
requires that the ultimate power source is star formation). 
Therefore the gas fuelling the burst is unlikely to be dis- 
tributed over the whole galaxy. Scenarios do exist in which 
positive feedback can happen where explosions in one part of 
a galaxy can trigger those in another part (e.g. Taniguchi, 
Trentham & Shioya 1998), but these mechanisms are not 
powerful enough to generate the kinds of luminosity required 
here over a whole galaxy. 

The important point is that if Al is correct, then the 
above suspicion of the connection between local ULIGs and 
the SCUBA sources based on the (remnant) number density 
and SED concordance now becomes a very strong assertion. 
This is because the only things that the SCUBA sources can 
evolve into are elliptical galaxy cores and their associated 
black holes. 

We do not consider the possibility that stars are made 
in the nucleus of the galaxies (as would be required by Al) 
and then redistributed throughout the galaxy in stellar- 
kinematical mergers because (1) the dense stellar cores so 
produced will not be disrupted by mergers (or secular evo- 
lution), (2) we would require about 100 mergers per L* ellip- 
tical and the observed merger rate (e.g. Carlberg et al. 2000) 
is not that high, and (3) the number density concordance de- 
scribed above is then lost due to the large number of merg- 
ers. 

Van der Marel 1999 finds the black hole mass (in so- 
lar units) log 10 Mbh ~ — 1-83 + log 10 Lv ■ Therefore for a 
Schechter (1976) L* elliptical galaxy (M v = -21.5) at red- 
shift 2 = 0, the central black hole has mass 5.1 x 1O 8 M . 
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For the same elliptical galaxy, the stellar core mass is 
M* c = 1.6 x 10 9 Mq, assuming that the core is approxi- 
mately isothermal at small radii (r < r c where r c is the core 
radius), that the galaxy has a de Vaucouleurs r 1 ^ 4 profile at 
large radii, that the ratio of r c to the de Vaucouleurs effective 
radius r e is 0.033 (derived from the sample of Lauer 1985), 
and a value of (M/Ly) = 3.86 for elliptical galaxies (see 
the discussion in Section 2.1.3 of Fukugita et al. 1998 - this 
assumes the stellar initial mass function (IMF) of Gould, 
Bahcall & Flynn 1996, which turns over at about 0.5 Mq, 
similar to the IMF of Kroupa, Tout & Gilmore (1993) that 
we adopt below). Therefore, for typical luminous elliptical 
galaxies 



(1) 



M BH _ 1 

M* c ~ 3 

Given the scatter in the correlations presented by van der 
Marel and Lauer, and the fact that many cores deviate sig- 
nificantly from isothermality, the most extreme galaxies may 
deviate from this ratio by as much an order of magnitude. 
For accretion onto a black hole, 



= r/Mc 



(2) 



where r\ is the efficiency of accretion and M is the mass 
accretion rate. 

For dust-enshrouded star formation where almost all 
the flux is absorbed and reradiated at far-infrared wave- 
lengths, the star formation rate (in Mq yr _1 ) is (Rowan- 
Robinson et al. 1997) 



1.2 x 10" 1U L* 



(3) 



where L, (in Lq) is the total luminosity generated by the 
young stars. This relation adopts a value of <p = 0.45 ap- 
propriate to the stellar IMF of Kroupa et al. (1993), and a 
value of e = 1 (see Section 3 of Rowan- Robinson et al. 1997 
for a definition of those parameters and a discussion of the 
derivation of this equation and other approaches) . The IMF 
chosen here has been shown to be consistent with most cur- 
rent observations of both star-forming regions and the solar 
neighbourhood (Gilmore & Howell 1998). 

Combining the two previous equations yields 



Lacc/L, ~ 1700 Tj. 



(4) 



Accretion is a lot more efficient than star-formation at gen- 
erating luminosity per unit mass: a 10 12 Lq galaxy has a 
star formation rate of about 80 Mq yr" 1 but a 2 x 10 1 Lq 
quasar has an accretion rate of about 1 M yr _1 for r) = 0.1 
(e.g. Chokshi & Turner 1992). Therefore averaged over time 
in the SCUBA sources: 



< L /L, > t ~ 560j? 



(5) 



This ratio is very big, meaning that the total power gen- 
erated in the SCUBA sources comes almost entirely from 
accretion if Al is correct. Note that throughout this calcu- 
lation, we have been working backwards, deriving the prop- 
erties of the progenitor from the properties of the remnants 
under the assumption that we can match particular progen- 
itors to particular remnants; this is quite a different calcula- 
tion from one in which we start out with a cloud of gas and 
attempt to track its progress. 

The fact that the power sources are so heavily accretion- 
dominated means that most SCUBA sources must be dust- 



enshrouded active galactic nuclei (AGNs) if Al is correct. 
This is not surprising, since the SCUBA sources contribute 
at least 10 per cent of the bolometric luminosity density of 
the Universe, and this could not be generated from the stars 
in gE cores alone, which comprise less than 1 per cent by 
mass of the stars in the Universe, and these are the only 
stars that the SCUBA sources are allowed to make if Al is 
correct. 



3 ADDITIONAL EVIDENCE FOR THE SCUBA 
SOURCE - AGN CONNECTION 

There exist other reasons for suggesting a connection be- 
tween the SCUBA sources and dust-enshrouded AGNs. 

1. The local supermassive BH density (Magorrian et al. 1998, 
van der Marel 1999) is high and its production generates a 
bolometric background of 



5.1 



PAGN, dusty 



2 x lO 5 M Mpc" 



nWm 2 sr 1 



(Trentham & Blain, 2000) if they radiate at one-tenth of 
the Eddington luminosity and if the luminosity density of 
obscured AGNs follows that measured by Boyle & Terlevich 
(1998) for optical quasars. This equation follows from a con- 
sideration of the energy released by accretion processes per 
comoving volume element given a final MDO density some 
part pagn, dusty of which was generated by dust-enshrouded 
accretion. The fiducial Pagn, dusty is the estimate of Salucci 
et al. (1999) for obscured quasars. The SCUBA sources in 
Fig. 1 generate about 7 nWm _2 sr _1 , which is close to this 
value. 

2. There is growing evidence that the hard (30 keV) X-ray 
background originates from absorbed quasars which reradi- 
ate energy absorbed at optical, ultraviolet, and soft X-ray 
at far-infrared and submillimetre wavelengths. Were this the 
case, Fabian & Iwasawa (1999) find a total reradiated energy 
density of about 3nWm~ 2 sr _1 . This is close to the numbers 
in the previous paragraph, and it is therefore plausible that 
these same AGNs that we are hypothesizing are the SCUBA 
sources are also the objects which contribute to the hard X- 
ray background (see also Almaini et al. 1999 and Bautz et 
al. 2000) . Note that if Compton-thick sources are very com- 
mon (see Maiolino et al. 1998), the reradiated energy will 
be somewhat higher than 3nWm _2 sr _1 . 

3. The Madau plot - that is, the cosmic star-formation rate 
of the Universe as a function of redshift - when determined 
by ultraviolet and optical measurements alone (e.g. Steidel 
et al. 1999) and integrated over time, produces a total stellar 
density in critical units of fi* ~ 0.004 (Pettini 1999), equal 
to the observed local value (Fukugita et al. 1998) Therefore 
an additional population of objects which are star-forming 
would lead to f2* being overproduced. Since > 10 per cent 
of the bolometric luminosity of the Universe is generated 
in the SCUBA sources, were they to be star-forming, they 
would be producing a substantial fraction of fi* . We there- 
fore require an alternative power source, such as accretion 
onto a black hole. The Madau plot described above includes 
a correction for star-formation obscured by dust in normal 
galaxies, which are not the same as the ULIGs or SCUBA 
sources described in the previous sections: their predicted 
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850-pm fluxes, even for the most luminous normal star- 
forming galaxies, are substantially lower (see e.g. Trentham 
& Blain 2000). 

4. All supermassive BHs at low z are found in spheroids 
(Magorrian et al. 1998), except those in ULIGs like Mrk 
231 which are probably in the process of forming spheroids 
(the core in the dissipative collapse, the rest of the galaxy 
in the stellar-kinematical merger of pre-existing stars). No 
supermassive BHs are known in pure disks. Such a scenario 
follows from the physical processes described in Section 2. 

5. In their spectroscopic survey, Barger et al. (1999) find 
some dust-enshrouded AGNs. But it is clear that powerful 
Seyfert 1 galaxies are not present in substantial number, so 
that if this sample is AGN dominated, they must be very 
highly obscured. Such objects are not very common in local 
ULIG samples (see Section 3.2 of Lutz, Veilleux & Genzel 
1999). 



4 PROBLEMS WITH THE SCUBA SOURCE - 
AGN CONNECTION 

The previous two sections have provided some evidence in 
support of a picture where Al is true. But there are prob- 
lems with such a scenario: 

1. The ratio of the number of ULIGs powered by accretion 
onto a BH to the number of ULIGs powered by star for- 
mation must change from about 1/3 at z = (Genzel et 
al. 1998) to about 80 at the redshift of the SCUBA sources 
(for r\ = 0.1) if Al is correct. This is very large effect, and 
there exists no clear physical mechanism to account for such 
a big difference, particularly since the low-redshift ULIGs 
and high-redshift SCUBA galaxies are similar in every other 
respect if Al is true. 

2. The q value (relating the far-infrared flux to the 4.85 
GHz radio flux; Condon, Frayer & Broderick 1991) of SMM 
J02399— 0136 is high, suggesting that about half of its power 
(at most 3/4) comes from star formation (Frayer et al. 1998). 
Recent X-ray measurements of this galaxy with Chandra 
(Bautz et al. 2000) also suggest that the fraction of the bolo- 
metric luminosity contributed by star-formation is at least 
20 %, probably more. This is much higher than the value of 
1/(560j7 + 1) predicted in Section 2. For this not to be a con- 
cern, SMM J02399— 0136 must be a very unrepresentative 
example of a SCUBA galaxy (which is not the case given 
the optical spectroscopy of Ivison et al. 1998 and Barger et 
al. 1999) or its AGN must emit an amount of radio lumi- 
nosity given its bolometric luminosity that somehow mimics 
starbursts. 

3. If the SCUBA sources are powered by AGNs, the infalling 
gas must lose at least 99 % of its initial angular momentum 
before reaching the BH event horizon. This is not achive- 
able by any known mechanism. Suggested mechanisms in- 
clude global gravitational or magnetic processes (Begelman 
1994), but the details are unclear. Note that the radius of the 
orbit of the two nuclear gas disks in Arp 220 is about 250 pc 
(Sakamoto et al. 1999), which is approaching the size scale 
on which the central black hole dominates the gravitational 
potential. Although this has been a long-standing problem, 
the increased angular resolution of the gas achieved by cur- 



rent inteferometry gives it an additional importance if Al 
is to be regarded as a plausible hypothesis. 

One possible solution that simultaneously addresses all 
three problems above could be that: the ULIGs undergo an 
intense star formation phase in the center that generates 
a bolometric luminosity of 10 12 Lq followed immediately by 
an AGN phase that lasts for a shorter time, but generates 
a higher bolometric luminosity, say about 10 13 Lq. If one 
selects objects instantaneously with bolometric luminosities 
above 10 12 Lq, one would end up with a sample dominated 
by starbursts, perhaps in a ratio similar to what Genzel et 
al. (1998) found locally if the relative timescales and lumi- 
nosities between the starburst and AGN phase are extremely 
large and small, respectively. This could remain true even if 
the time-averaged luminosity was generated predominantly 
by AGNs, as required by eq. (5). Such a scenario could help 
the angular momentum problem since the black hole could 
accrete the stars themselves (Begelman 1994); given eq. (1) 
about one-fourth by mass of the total stars produced would 
have to be accreted. It is also consistent with a picture in 
which the most luminous things that exist at any redshift 
are AGN-powered (see e.g. Sanders & Mirabel 1996). 

Another possible solution is that the fraction of ULIGs 
at z — that are AGNs has been underestimated because 
the AGNs are so heavily obscured that they are invisible 
even at mid-infrared wavelengths. The detection of hard X- 
ray emission from NGC 6240 (Vignati et al. 1999) suggests 
that the bolometric luminosity of this infrared-luminous 
galaxy, previously thought to be starburst-powered, orig- 
inates from an AGN. If this object is at all typical of 
ULIGs, then the possibility of a much higher ratio of local 
ULIGs that are AGN-powered than one-third seems possi- 
ble (see however the discussion in the final section of Lutz 
et al. 1999). 



5 DISCUSSION AND SUMMARY 

We have shown that if the SCUBA sources, which appear 
on SED and luminosity grounds to be ULIGs at high z have 
the same gas density structure as local ULIGs, then they 
must generate most of their power from embedded AGNs. 
Such a scenario is attractive for many reasons: it generates 
the correct local number density of elliptical galaxy cores 
of the correct mass, the correct local mass density of their 
associated black holes, and gives a formation mechanism 
for kinematically decoupled cores (which probably form in 
some other way from the rest of the galaxy). These SCUBA 
sources would then be very different from local ULIGs, which 
derive most of their power from star formation. 

But there are problems (see Section 4) . If these cannot 
be resolved, then Al must be abandoned. In this event, the 
gas in the SCUBA sources must be far more diffuse than 
in local ULIGs because in order to generate such huge lu- 
minosities, a substantial fraction of the stars in the galaxies 
must be made; in z = galaxies, large fractions of the stars 
in galaxies do not reside at extremely high densities of 100 
Mq pc -3 or greater so it is improbable that they were made 
out of gas at these densities. In local ULIGs all the bolomet- 
ric luminosity seems to come out of regions associated with 
gas at this high density. 

Either way, there is a big difference between local ULIGs 
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and high-redshift SCUBA sources. This is the main result 
of this letter. 

Testing the validity of Al will be an important exercise 
in the future. HCN measurements (like for the Cloverleaf 
quasar) would be best, but until these and other millimetre 
wave spectral diagonstics (other than CO) are available, the 
best hope is to rely on indirect methods like those outlined 
in the letter. 
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